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1.  INTRODUCTION 


An  important  aspect  of  the  electromagnetic  pulse  (EMP)  vulnerability  assessment  of 
electronic  equipment  is  the  question  of  whether  a particular  semiconductor  component  will  be 
damaged  by  the  EMP-induced  transient.  The  methods  currently  used  to  predict  semiconductor 
damage  are  based  on  calculated  power  dissipation  in  the  device  junctions.  Most  such  efforts 
involve  replacing  the  highest  peak  of  the  power-versus-time  curve  by  a rectangular  pulse  which 
is  approximately  equivalent  to  it  in  energy  content.  Then  the  damage  data  for  that  device  are 
scanned  to  see  if  the  calculated  power-width  coordinates  for  the  pulse  fall  above  or  below  the 
nominal  damage  curve.  Vulnerability  assessments  in  such  a scheme  are  quite  imprecise,  usually 
reducing  to  such  statements  as  the  system  is  “vulnerable”  or  is  “not  vulnerable.” 

This  procedure  appears  reasonable  enough,  because  any  such  analysis  must  be  based  on  the 
existing  device  damage  data,  and  these  were  almost  invariably  obtained  by  step-stressing  a 
device  to  failure,  using  rectangular  pulses.  It  has,  however,  two  serious  faults;  (1)  the  treatment 
of  arbitrary  power-versus-time  waveshapes  as  rectangular  pulses  and  (2)  the  lack  of  a probabilis- 
tic treatment. 

This  report  describes  two  possible  improvements  on  the  above  procedure.  The  first 
improvement  addresses  the  statistical  distribution  of  damage  by  assuming  that  the  device  damage 
power  is  lognormally  distributed  at  constant  pulse  width.  The  nominal  damage  level  is  then  used 
with  an  estimate  of  the  standard  deviation  in  power  to  damage,  after  review  of  a large  amount  of 
data,'  to  calculate  the  probability  of  device  failure,  given  power  inputs  in  the  form  of  rectangular 
pulses.  The  other  improvement,  which  removes  both  of  the  above  objections,  uses  the  entire 
power-versus-time  curve  directly  to  calculate  the  probability  of  device  failure. 

2.  DEVICE  DAMAGE  DATA  DISTRIBUTION 

Device  damage  data  are  generally  obtained  by  step-stressing  a number  of  devices  to  failure, 
using  rectangular  pulses  of  various  widths.  The  results  of  this  procedure,  after  proper  curve 
fitting,  provide  the  device  damage  curve — a log-log  plot  of  nominal  damage  power,  P^,  versus 
pulse  width,  t.  There  is  substantial  disagreement  about  the  form  of  Po,  mostly  due  to  the  lack  of 
good  experimental  data,  and  it  is  likely  that  no  single  form  can  adequately  describe  all  devices. 
Still,  for  very  short  pulses,  the  adiabatic  approximation*’®  appears  generally  valid.  That  is 

Poir)  « T-', 

while  for  very  long  pulses  thermal  equilibrium  can  be  reached,  that  is. 

Pair)  a T*. 

For  the  intermediate  region,  there  is  a general  tendency  to  use  the  Wunsch-BelP  model,  with 

Po(r)  “ r-*'*. 

It  is  not  at  all  clear  where  these  regions  should  meet.  The  available  data  imply  that  the 
adiabatic  approximation  is  excellent,  at  least  up  to  50  to  100  ns,  while  the  equilibrium 


' Defense  Nuclear  Agency  EMP  Handbook,  Component  EMP  Sensitivity  and  System  Upset,  Ch.  13  (22  September 
1975). 

' D.  C.  Wunsch  and  L.  Marzitelli,  Semiconductor  and  Nonsemiconductor  Damage  Study,  Braddock,  Dunn,  and 
MacDonald,  Inc.  (April  1969). 

* D.  M.  Tasca,  IEEE  Trans.  Nucl.  Sci.,  NS-17,  6 (December  1970). 

* D.  C.  Wunsch  and  R.  R.  Bell,  IEEE  Trans.  Nucl.  Sci.,  NS-15,  6 (December  1968). 


approximation  is  usable  down  to  perhaps  a few  seconds.  Between  these  limits  is  a broad 
transition  region  wherein  Po  varies  smoothly  from  one  limit  to  the  other.  On  log-log  paper,  the 
slope  of  Pd  varies  monotonically  from  - 1 (adiabatic)  to  0 (equilibrium).  The  Wunsch-Bell  model 
represents  the  transition  region  by  a line  of  slope  -Vi. 

Most  cases  of  interest  to  us  involve  total  times  less  than  about  500  ns.  Thus,  the  equilibrium 
region,  and  most  of  the  transition  region,  as  well,  can  generally  be  ignored.  For  this  reason,  and 
because  of  the  large  amount  of  scatter  in  the  damage  data,  it  appears  practical  to  approximate  Pd 
by  two  straight  lines  which  intersect  at  a certain  point,  namely. 


_k'T-‘,  T<T, 
fD  - .r  a:  To 


where  K is  the  Wunsch-Bell  damage  constant,  and  K'  is  determined  by  the  requirement  that  the 
two  lines  intersect  at  t = To . That  is, 

/’z)(to)  = A’Vo"*  = ATto'"* 


K'  = jICto''*. 

In  a later  section,  we  shall  represent  the  probability  of  device  damage,  for  a power  input 
pulse  of  amplitude,  W,  as  a cumulative  lognormal  distribution: 


^0(1^)  = Vi'  1 + erf 


InW-M 

<r\/2 


where  = InPo,  and  (r  is  the  standard  deviation  in  the  natural  log  of  the  damage  power 
distribution. 

Examination  of  many  sets  of  damage  data*  has  shown  that  the  distribution  in  power  to 
damage,  at  fixed  pulse  width,  is  approximately  lognormal  and  that  the  typical  spread  is  about  a 
factor  of  five  in  power.  Thus,  in  the  lognormal  description^  the  standard  deviation,  (t,  is 
approximated  by 

<r  = ln5  = 1.61  (see  app  B). 

The  median  value  of  damage  power  is  given  by  Pd-  Therefore,  the  mean,  n,  of  the  damage 
distribution  is 


_ /ln(A:'T  *),  adi 
^ lln(^T-‘'*),  Wu 


adiabatic 

Wunsch-Bell 


These  two  moments  completely  describe  the  distribution.  In  practice,  since  damage  data  are 
seldom  available  for  pulse  widths  less  than  50  ns,  we  find  K by  finding  the  line  of  slope  -Vi 
which  best  fits  the  data.  The  value  of  K'  is  then  determined  by  the  choice  of  Tg.  In  general,  we 
used  Tg  = 50  ns. 

3.  THERMAL  RESPONSE 

It  is  well  known  that,  if  one  obtains  the  response  of  a system  to  a step  excitation,  the 
impulse  response  can  be  gotten  by  differentiation,  and  that  the  response  to  any  other  excitation 
can  be  calculated  by  convolving  the  impulse  response  with  the  arbitrary  excitation.  Thus, 


' Defense  Nuclear  Agency  EMP  Handbook,  Component  EMP  Sensitivity  and  System  Upset.  Ch.  13  (22  September 


consider  an  initially  relaxed  system  described  in  terms  of  an  impulse  excitation,  6(t  - /')>  and  a 
response,  h(t  - /'): 

H-\t)h(t  - t')  = S(f  - t'). 

Here,  is  a linear  differential  operator  with  constant  coefficients.  If  we  assume  that  no 

signals  are  applied  before  f = 0,  this  expression  can  be  integrated  once  to  give 


f hit  - n dt'  = f 8(t-  n dt'. 

•'ft  •'ft 


Since  no  response  is  possible  before  an  excitation  is  applied,  the  integration  can  be 
terminated  at  t'  = t.  Thus, 

H-\t)  f hit  - t')  dt'  = f 8(f  - /')  df  = Uit), 

•'o  •'0 

where  Uit)  is  the  unit  step  function  defined  by 

Uit)A 


The  integral  on  the  left  now  gives  the  time-domain  response  to  the  unit  step.  Laplace 
transforming  both  sides  gives 


and 


r hit  - 1')  dt' 
f hit  - t')  dt' 

•'ft 


His) 

s 


= at^.(f)]  = v„(j). 


Then  His)  = sV;,(r),  where  His)  is  the  transfer  function  and  V«(j)  is  the  step-function  response  in 
the  frequency  domain. 


by 


The  impulse  response,  which  is  the  inverse  transform  of  the  transfer  function,  is  then  given 

hit)  = L-'[His)]  = 


Mf)  = ^Mf)]  + M0+)6W. 

The  last  equation  follows  from  the  two  standard  transforms: 


where 


L[f^\  = sFis)  -fi0+) 

^8(0]  = 1. 


Fis)  = Of  it)]. 
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The  response  to  a general  excitation,  e(t),  can  now  be  written 


v(t)=  f h(t  - t')e(t')  dt' 

= jf‘  { 5(737)  vjo+)s(t  - /')}e(f')  dt’ 

v(t)  = D,(0+)e(/)  + f e(t')  ~ - r')]  dt'. 

Jo  d{t  - t) 

In  the  application  to  thermal  damage  to  semiconductors,  we  make  the  following  equiva- 
lences. 

v{t)  -*  ^T{t)  = temperature  rise 
e{t)  -*  P{t)  = power  dissipation 
v^{t)  -*  R(t)  = response  to  rectangular  pulse  power 

ATW  = /?(0+)P(/)  + r Pin  ..  R{t  - n dt'. 

J0  a{i  t ) 

The  first  term  can  be  ignored  on  physical  grounds.  That  is,  the  response  at  zero  time  to  an 
excitation  which  began  at  zero  time  can  be  nonzero  only  for  inflnitely  large  power  input.  Hence, 

mt)  = I'  Pit') 

Now,  let  us  assume  tnat  the  temperature  rise  at  the  semiconductor  junction — due  to 
application  of  a rectangular  pulse  of  unit  amplitude  and  width  r — is 

ATfr)  = R(t). 

For  a linear  system,  the  temperature  rise  due  to  any  other  rectangular  pulse  of  the  same  width  is 
proportional  to  Rir).  In  particular,  for  a pulse  equal  to  the  nominal  damage  power  for  that  width, 
we  expect 

ATXt)  = Tc=  Poir)Rir), 

where  Tc  is  the  temperature  at  which  damage  occurs.  Thus,  we  can  write 

Dividing  by  Tc  gives 


This  is  a very  convenient  relation  between  the  instantaneous  junction  power  dissipation 
and  the  normalized  temperature  rise  in  the  junction.  As  long  as  / < 1,  we  expect  no  damage. 
When  / =■  1,  we  usually  make  the  nominal  prediction  of  damage.  In  the  next  section,  we  shall  see 
how  these  predictions  can  be  made  more  quantitative. 


Substituting  the  two  straight-line  expressions  for  Pj, , which  we  exhibited  in  section  2,  gives 


f‘‘  d (t-  t')*'* 


/,  2a:(i ' '' 


The  function  h(t)  described  earlier  in  this  section  will  be  recognized  as  the  Green’s  function 
of  the  operator  H“‘(t). 

4.  DAMAGE  ASSESSMENT 

We  now  have  in  hand  the  necessary  mathematical  apparatus  for  damage  assessment. 
Basically,  all  that  remains  is  to  exhibit  the  answers  for  the  two  methods  mentioned  in  section  1. 

In  the  first  method,  we  assume  power  inputs  in  rectangular  pulse  form.  Since  we  know  that 
the  device  damage  data  are  lognormally  distributed,  we  can  immediately  write  down  the 
probability  that  a pulse  of  amplitude  W and  width  r will  cause  failure.  It  is  given  by  the 
distribution  function,  F.  Thus, 


where  n = InFo  and  o-  = 1.61. 

A simple  change  of  variable  gives 


F„(W)  = " 

v2ff 

Fo(W)  = Vi  1 -t-  erf|- 


InW-fi 

0V2 


The  survival  probability  is  then  given  by 


P,{W)  = 1 - FoiW)  = 1 


In  the  second  method  we  are  concerned  with  the  distribution  in  the  maximum  values  of  the 
convolution  integral  7(0,  in  those  cases  where  damage  results.  For  rectangular  pulses,  it  is 
obvious  that  the  distribution  in  this  /„„  is  lognormal  with  the  same  a-  as  the  pulse  power.  If  the 
thermal  model  is  to  be  useful,  this  must  also  be  true  with  arbitrary  power  waveforms,  because 
the  distribution  is  really  in  junction  temperature.  Thus,  since  the  median  value  of  /„ax  which 
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or 


P,  = F(0.181)  = 0.572. 


5.2  Rectangular  Pulse  Case 
5.2.1  Adiabatic  Model 

First,  we  compute  the  most  probable  value  of  damage  power  for  pulse  width  r,  using  the 
formula  appropriate  for  the  adiabatic  model. 

M = ln{^’T-')  = ln(A:T„>'*T->) 

M = In^  + VilnTj  - Inr 

fji  = ln0.91  + Viln5  + 41nl0  - lnl.36 

H = -0.09  + 0.80  - 0.31  + 9.21 

/X  = 9.61. 

In  addition,  we  have 

InW  = ln9840  = 9.19. 

• There  is  a slight  change  in  notation  here.  What  we  previously  called  Fo(/)  is  really  This  distinction  is 

necessary  for  use  of  probability  tables. 
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Putting  these  into  the  formula  for  survival  probability  gives 

or 

P,  = 0.60. 

5.2.2  Wunsch-Bdl  Model 

This  calculation  is  entirely  analogous  to  the  adiabatic  calculation.  Thus, 

M = = InA:  - Vilnr 

fi  = ln0.91  - Vilnl.36  + 41nl0 

At  = 8.81, 

and  the  survival  probability  is 

or 

P,  =0.41. 

For  this  particular  problem,  the  convolutional  approach  gives  a survival  probability  which  is 
about  5 percent  lower  than  that  given  by  the  adiabatic  model,  but  about  39  percent  higher  than 
that  given  by  the  Wunsch-Bell  model.  This  is  not  necessarily  true  in  other  cases.  It  is  easy  to 
construct  waveforms  having  multiple  peaks  close  ^together,  for  which  the  convolution  integral — 
since  it  takes  account  of  all  peaks — would  give  a significantly  lower  probability  of  survival  than 
would  either  of  the  other  two  models. 


6.  DAMTRAC  MODIFICATIONS 

One  major  component  of  the  machinery  necessary  for  damage  assessment  is  a computer 
program  which  calculates  the  response  of  a circuit  to  input  transients.  We  used  an  extended 
version  of  a program  called  DAMTRAC,*  a version  of  the  TRAC  code,  which  calculates  power 
dissipation  in  given  circuit  components  as  a function  of  time. 

Almost  all  of  DAMTRAC  is  kept  in  object  form  on  a permanent  file.  Subroutine  TRAEQ  is 
available,  however,  and  this  is  where  the  necessary  changes  were  inserted.  Subroutine  TRAEQ 
is  called  several  times  for  each  solution  of  the  circuit  equations,  and  the  current  values  of 
elapsed  time  and  power  dissipation  are  then  available.  Each  time  an  acceptable  solution  of  the 
circuit  equations  is  reached,  the  current  values  of  the  convolution  integrals  are  calculated.  At  the 
end  of  a predetermined  interval,  the  maximum  values  of  these  integrals  are  found  and  the 
probabilities  of  component  failure  are  calculated.  At  the  end  of  the  job,  the  integrals  are  plotted 
as  functions  of  time  by  specifying  appropriate  node  numbers  in  the  input  data  deck.  Up  to  10 
semiconductor  components  can  be  analyzed  for  each  DAMTRAC  run.  If  necessary,  the  program 
can  easily  be  modified  to  accommodate  a larger  number  of  devices.  A listing  of  this  subroutine  is 
given  in  appendix  A. 

*G.  H.  Baker,  A.  D.  McNutt,  D.  M.  Rubenstein,  and  G.  B.  Shea,  Damage  Analysis  Modified  TRAC  Computer 
Program  (DAMTRAC),  Harry  Diamond  Laboratories  TM-75-6  (May  1975). 
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7.  CONCLUSIONS 

Two  improvements  in  EMP  circuit  damage-analysis  techniques  have  been  described:  (1)  use 
of  a convolution  integral  to  account  for  arbitrary  power  pulse  shapes,  and  (2)  use  of  a routine  to 
give  probability  of  damage  rather  than  purely  qualitative  results.  These  changes  result  in  a 
convenient  calculational  tool  for  damage  assessment.  The  increase  in  running  time  over  the 
present  DAMTRAC,  with  the  rectangular  pulse  approach,  is  small.  This  method  is  substantially 
more  satisfying  because  of  the  quantitative  damage  probability  results  and  because  the  entire 
waveform  is  taken  into  account.  This  method  will  be  used  as  Harry  Diamond  Laboratories 
primary  circuit-analysis  technique  for  the  INCA  program. 


APPENDIX  A. 

SUBROUTINE  TRAEQ  LISTING 


A listing  of  the  subroutine  TRAEQ  follows.  TRAEQ  is  called  several  times  for  each  solution 
of  the  circuit  equations,  and  the  current  values  of  elapsed  time  and  power  dissipation  are  then 
available. 


C SUBROUTINE  TRAEQ 

subroutine  TRAEQ(KKA) 

implicit  HEAL*b  <A-HtO-/:)t  IVTEUER*2  (I-N) 

REAL*4  ST.VP»XP,ABC.E04«Tf 
InTESER*2  err 
InTEGER*4  KK4 

INTE6£R*4  CONPtCONO*CONT*LIMlT*KSTRT 
DIMENSION  DAMP ( 1 000 t 10) t P4R (1 000 • 1 0 ) > T IME ( 1 000 ) 

DIMENSION  CO (10) tUAMMAX(lO) «P0AM( 10) 

COMMON/ALPHA/A(100) •E(lVt40t2)«H(100tl00),O«T(lO0) tV (e40) « AL (SO) t 
IBVf  CT(11«2)«  01(30)t  DM«  0R(30)«  DV(30)t  El(20)«  E2(19)« 
2TEtTl«TOtVS(20t5) «V1 ( 240 ) • V2 ( 240 ) ,UC1 (30) « 

38CR(30) tBCV(30) tBEI (30) •c)ER(30) tBE V ( 30 ) «DIP ( 30 ) tEPl ( 20 ) «EP2 (20 ) t 
4UCR(30*7) tSlltTEXtTEltTOOdO) t TUP t TOl t VER ( 9) tAAAA(60) tBCDTt 
SBC  IP (30) tBCUR(290) tBE IP (30) tOTlMtORSPt 

7PpIC(30)t  PPIO(30)t  PPlE(30)t  PMCR(30)f  Pi«TR(30)t  OTAN(30tlb)t 
BRMAGt  TEllt  TOLlt  T0L2t  ruT2t  TUllt  VALl(200)t  VAL2(200)t 
9VAL3(200)t  BCUTlt  BCURl(290)t  PPlCl(30)t  PPIDl(30)t  PPIEl(30)t 
ASYMBltSYMB2tTITLL(9) tTOLPLtT0MINiST(19t2) tABCt 
•E04tCONO(30tS) tCONP(200tB) tC0NT(30tS) tERRt lAMi JSJt lOCt NCUTt 
BKTOTt  IBCt  NAMt  ITOTi  KTUTlt  NVt  Nlit  JJt  N3Bt  NBRANt  NBLOCKt 
CITOTlt  JUJt  N2Bt  ISTRKt  lSYM4t  NZt  NUt  NBBt  NvMt  NCRt  NNNNlt  ICPt 
OIPRINTt  NTPLPt  NPP(60)t  NOP(4)t  NAt  NZl<  NOlt  JJJt  KEi  NTRt 
EIFGOTCt  JVt  IPXt  NAlt  NTt  NBt  NlBt  JCATt  NCORtt  UVJt  NPt  IPLCONt 
FISYH3t  NPl(lUt2)t  J6(19) 

COMMON/8ETA/VP(5000) fXP(700) tTF (3b0) 

COMMON/GAMMA/PmCRFO(30) tPwCRRV(30) tBRKVOL(30) tOAMNK(30) 
COMMON/OELTA/PMTEFO(30)  tP«(TERV(30)  tP*)TCF0(30)  tPwTCRV(30) 

DATA  M/l/t  NPR/0/ 

C FOR  SUPPRESION  OF  M*T  MATRIX  PRINT  , SET  ERR=2 
ERR*2 

GO  TO  (9000t9003t9002t90Ul ) t KKA 

9000  CONTINUE 

C PART  1 - MATRIX  AND  TIMt  FUNCTION  EU''S.  AFTER  THIS  CArD. 

El (3)«V(21)*.42 

OMEGArZ. *3. 1416*30. OOE^Ob 

El (2) »9.*USIN(0ME6A*TE) ♦!.£-& 

C RETURN 

9003  CONTINUE 
C PART  2 

Rh  TURN 

9002  CONTINUE 

return 

9001  CONTINUE 

C part  4 - AUXILIARY  EQUATIONS  AFTER  THIS  CARO. 

V(24)sPhCRF0(1) 

V(2S)«PWCRRV(1) 

V(2b)>PMTCFO(l) 

V(27)*PMTCRV(1) 

V(2«)*P«*TEFD(1) 

V(29)>PMTERV(1) 

V(30)bPWCRFD(2) 

V(31 )bPMCRRV(2) 
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o o n n r>  n n 


the  FULL0MIN6  CAHOS  ARE  FOR  OARAGE  CALCULATIONS 

CO  - MUNSCH  DAMAGE  CONSTANT 
TO  • CROSSOVER  WlOTHi  hUNSCN  TU  AOIAOATIC  MODEL 
NDEV  - NUMBER  OF  VULNERABLE  DEVICES 

T0«5.E-8 

N0EV»2 

CO(I)>0.91 

CO(2>*U.20 

IF(TE.LE.O..OR.M.GE.1000)  SO  TO  200 
PMR<Ntl)>PWTCRV(l> 

P«R(M«2)>PHTERV(1) 

TIME<M)»TE 

IF(M.EO.l)  60  TO  190 

IF(T1MC(M) .LE.TIME(M-l) ) 60  Tu  200 

00  10  1>1*NDEV 
OAMP(M»I)>0. 

00  10  J»2*M 

IF(T1ML(M)-TIME<J)-T0)  Stbtb 

5 DAMP(M»I)»0AMP<M»I)* (TIME! J) -TIME (J-l» )*(PRR(J,1) ♦PBRCJ-ltl) )/2./ 

1 DSURT(T0)/CD(I) 

60  TO  9 

b OAMP(MtI)xOAMP(M«I)« (T1ME( J)-TIME(J-l) )*(PWR( J«I>/DSORT(TIME(M>- 
1 TIME(J) )*<P«(R(J-Itl)/DSURT(TIME(M)-TIME(J-1) > ) )/A./CD<I> 

9 V(100»I)xOAMP(MtI) 

10  CONTINUE 

IF(NPR.EQ.I)  60  TO  190 
lF(TE.6E.CT(lt2r-CT(lfl) ) SO  TO  ISO 
60  TO  190 
ISO  NPRxl 

170  DO  IBO  IxltNOEV 

0AMMAX(I)30AMP(2tl) 

00  17s  J«3fM 

17S  lF<OAMP(JtI)  .OT.OAMMAXd)  ) OAMMAX  (1 ) sOAMP  ( I ) 

Ys0L06(0AMMAX(l) )/l.bl/SURT(2.) 

IBO  POAMII)x0.5*(l.«OERF(y)) 

HRITE(b«16S) 

IBS  FORMAT(1H1«SOX«30HSUMMARY  OF  DAMAGE  CALCULATIONS/) 

WRlTE(b*IBb)  <I«OAMMAX(I) tlxl.NOEV) 

iBb  FORMAT (//25X*4lHMAXIMUM  VALUE  OF  CONVOLUTION  INTEGRAL  N0.tI3t2H  «f 
1 1PE11.3) 

NRlTE(btI87)  (ItPDAM<I) «I>1«N0EV) 

187  F0RMAT(//2SX«35HPR0BAblLITY  OF  DAMAGE  TO  DEVICE  N0.tI3«2H 
1 IPE11.3) 

MRlTE(btI88> 

188  FoRHATdHl) 

190  MsM*l 

200  return 

END 
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APPENDIX  B 

A NOTE  ON  SEMICONDUCTOR  DAMAGE  DATA 

The  value  assumed  for  <r  in  section  2,  body  of  report,  may  be  excessively  large  in  most 
cases.  If  so.  the  reader  is  free  to  substitute  any  value  he  prefers.  It  should  be  pointed  out. 
however,  that  semiconductor  damage  tests,  in  general,  were  not  statistically  designed.  That  is, 
the  experimenter  made  no  effort  to  plan  his  tests  so  that  the  results  would  be  statistically 
meanin^l.  Sample  estimates  of  <r  and,  to  a lesser  extent,  n,  based  on  such  tests  may  be 
misleading  to  an  engineer  and  downright  useless  to  a statistician. 
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